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Abstract 


It  is  argued  that  only  phase,  modulo-2ir.  Is  of  Interest  In  data 
communication  applications.  The  probabilistic  evolution  of  random  walk 
phase,  modulo-2ir.  Is  studied  and  the  results  used  to  derive  a MAP  se- 
quence estimator  for  phase  on  the  Interval  [-iT,7r).  The  MAP  sequence 
estimator  Is  a Vlterbl  tracker  (or  dynamic  programming  algorithm)  for 
tracking  phase  on  a f Inlte-dlmenslonal  trellis  In  [-Tr,iT).  Computational 
requirements  are  discussed  and  several  representative  simulations 
are  Included.  Applications  of  the  algorithm  to  carrier  phase  tracking 
In  coherent  data  communication  systems  are  outlined. 


I.  Introduction 


A general  model  for  data  transmission  over  an  Imperfect  bandllmited 
channel  Includes  the  effects  of  Intersymbol  Interference,  carrier  phase 
fluctuation,  and  additive  receiver  noise.  The  carrier  phase  fluctuation 
arises  from  changes  In  the  channel  transmission  characteristics,  ripple 
In  dc  power  supplies,  and  timing  Jitter  at  the  transmitter  and  receiver. 

A somewhat  primitive,  but  tractable,  model  for  carrier  phase  fluctuation 
Is  the  so-called  random  walk.  A more  sophisticated  model  would  account 
for  the  fundamental  and  harmonic  components  of  ringing  currents  and 
60Hz  Interference  [1]. 

In  this  paper  we  propose  a new  approach  to  carrier  phase  tracking. 

We  observe  that  the  measured  signal  In  a data  communication  system  is 
Invariant  to  modulo-2Ti  transformations  on  the  carrier  phase.  Therefore 
maximum  a postevioH,  (MAP)  sequence  estimation  for  a modulo-2ir  version 
of  the  underlying  phase  sequence  seems  an  appropriate  estimation  problem 
to  pose.  To  solve  this  sequence  estimation  problem  we  derive  an  expres- 
sion for  transition  probabilities  In  a modulo-27r  version  of  random  walk. 

A rigorous  treatment  may  be  given  by  studying  random  walk  on  the  circle 
[2].  The  modulo-2ir  transition  probabilities  are  used  to  derive  a path 
metric  for  use  In  a Vlterbl  algorithm  for  tracking  phase,  modulo-2‘tr,  on 
a flnlte-dlmenslonal  trellis  In  [-■ir,ir).  The  resulting  algorithm  approxi- 
mates the  maximum  a posteriori  (MAP)  modulo-2ir  phase  sequence.  Several 
representative  simulations  are  presented  to  Illustrate  the  tracking 
characteristics  of  the  proposed  algorithm.  Applications  of  the  algorithm 
are  outlined  for  phase  coherent  data  coimnunlcatlon  over  channels  that 
may  or  may  not  have  been  equalized. 


m 
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II.  Signal  and  Phase  Models 

We  begin  our  development  by  considering  Che  output  of  a quadrature 
demodulator  when  there  Is  no  data  transmission  taking  place.  This  allows 
us  to  formulate  and  solve  the  modulo-2Tr  carrier  phase  tracking  problem. 

The  solution  Is  generalized  In  Section  Vlll  to  Include  simultaneous 
carrier  phase  tracking  and  symbol  detection. 

Consider  the  following  complex  representation  for  the  sample  values 
appearing  at  the  output  of  a complex  demodulator: 

- exp(ji|>j^)  + n^  , k»0,±l,±2,  ...  (1) 

The  additive  noise  sequence  Is  assumed  to  be  a sequence  of  inde- 

pendent Identically-distributed  (l.l.d.)  complex  normal  random  variables: 

Uj^  : N(0,a2)  v^^  : N(O.a^)  (2) 

“k 

2 

The  notation  x:N(p,a  ) Indicates  that  x is  a real-valued  normal  random 

2 

variable  with  mean  p and  variance  o . The  notation  x J_L  y means  x and  y 
are  statistically  Independent. 

The  random  phase  sequence  appearing  In  (1)  Is  assumed  to  be  a 

random  walk  sequence  characterized  by  l.l.d.  normal  Increments  Wj^: 

*k+l  = ♦k  "k  (3) 

w^  : N(0,a^)  11  Wj^  V k 1 

Let  f(x/y)  denote  the  conditional  density  function  for  the  random 
variable  x,  conditioned  on  (a  realization  of)  y.^  Then  the  conditional 


''What  we  really  mean  here  Is  the  conditional  density  function  for  a 
random  variable  x,  given  that  y equals  some  value  a.  A more  appropri- 
ate, but  also  more  cumbersome,  notation  would  be  f^^  (‘/a). 


3 


density  for  given  la 

By  this  notation  we  mean  "...  conditioned  on  the  value  for  the  phase 

at  time  k,  the  phase  at  time  k+1  is  a normal  random  variable  with  mean 
2 

*,  and  variance  o ." 
k w 

Note  that  the  measurement  model  of  (1)  is  Invariant  to  a modulo-Zn 
transformation  of  the  phase  That  Is 

Zj^  » exp(jij^)  + nj^  (5) 

where  Is  the  modulo-Zn  version  of 

+ 12-n  some  1 ; e[-n,ir)  (6) 

This  observation  motivates  us  to  model  the  probabilistic  evolution 
of  the  phase  sequence  and  derive  a MAP  sequence  estimator  for  the 

modulo-2ir  sequence 

III.  Review  of  Ungerboeck  [1]  and  Magee  [3] 

Ungerboeck  [1]  and  Magee  [3]  have  treated  carrier  phase  tracking 
In  the  context  of  simultaneous  phase  estimation  and  symbol  detection. 
Ungerboeck  assumed  a perfectly  equalized  channel  and  Magee  allowed  for 
a finite-length  channel  Impulse  response.  For  purposes  of  comparison 
with  the  approach  proposed  In  this  paper,  we  present  a brief  review 
of  [1]  and  [3]. 

In  his  treatment  of  phase  sequence  estimation  Ungerboeck  observed 
that  the  phase  Is  a continuous-valued  variable  and  that  no  Vlterbl- 
llke  algorithm  was  directly  applicable.  He,  therefore,  discretized  the 
range  space  for  to  a grid  of  spacing  (In  our  notation).  He  called 
this  approximating  sequence  a sequence  with  discrete  binary  Increments 
and  the  corresponding  approximation  a delta-mod  approximation.  The 


A 


underlying  parameter  space  Is  now  countably,  rather  than  noncountably , 
infinite.  Ungerboeck  handles  this  difficulty  by  restricting  attention 
to  only  those  phase  states  that  lie  in  a region  around  the  developing 
most-likely  path.  Magee  does  not  really  implement  a phase  sequence 
estimator.  Rather,  in  the  interest  of  conserving  computation,  he  first 
minimizes  a path  metric  with  respect  to  the  current  phase  value  (at 
each  step  in  the  recursion)  and  then  uses  a d3mamic  programming  algorithm 
to  find  a MAP  data  sequence  estimate. 


IV.  The  MAP  Sequence  Estimation  Problem  for  Modulo-2'ir  Phase 
Consider  the  following  maximization  problem  with  respect  to  the 

_ g 

modulo-Zir  phase  sequence  • 

max  (7) 

Here  f(*,*)  is  the  joint  density  function  for  the  K measurements 

and  the  K modulo-2ir  phase  values  Maximization  of  this  joint  density 

function  is  equivalent  to  maximization  of  the  a posteriori  density 


{Zk>^). 


Write  the  joint  density  of  (7)  as 


- K l^k  2 

IT  N (e  *,o^) 

^ ^ k-1  \ 


(8) 


The  last  simplification  in  (8)  follows  from  the  fact  that  the  n^  in  (1) 
are  i.i.d.  normal  random  variables  and  the  fact  that  conditionally. 


(e  ^ , ah 


""k  2 

We  have  used  the  subscripted  notation  N (e  , a ) to  indicate  that 

*k 
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the  conditional  density  of  the  complex  random  variable  z.  (conditioned 

2 

on  is  normal  with  mean  e and  variance  o : 

9 9—1  1 ■ 9 

N (e  , 0 ) = (2x0  ) exp{ =■  \z  -e  (9) 

20^  ^ 

Dropping  phase  Independent  terms  we  may  write  the  MAP  sequence 
estimation  problem  as 

- K 1 •^♦k 

maxdnf  ({<<)  2Re  Z z e (10) 

,r  ,K  ^ 20^  k-1 

Or,  defining  the  envelope  statistic  Cj^  = |zj^|  and  the  phase  statistic 
= arg  z^  (so  that  = |zj^|  exp  ( j ) , we  may  write 

max  («,nf  ({^j^}^)  + ^ ^ \ cos(ij)j^  - (11) 

/T  iK  2a  k=l 

• K 

The  function  density  function  of  the  K modulo-2x 

phase  values  • We  study  this  object  in  the  next  section.  We  empha- 

size that  “ arg  z^  takes  on  values  only  in  l-x,x). 

V.  Probabilistic  Evolution  of  the  Modulo-2x  Phase  Sequence 
We  are  Interested  in  the  probabilistic  evolution  of  the  sequence 
Accordingly,  we  seek  an  expression  for  the  conditional  density 

This  will  permit  us  to  write 

-kK__ 

f({*^^}p  • ir  f(VVl^  ’ (13) 

k*l 

because  the  modu]o-2x  transformation  does  not  destroy  the  Markov  property 
property  In  a random  walk. 

Note  ™ay  be  written 


•^k+l  ^k+1  ’ ’k+1 


K + Wl 
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The  conditional  density  of  Is 

f<Vi/V  ; N(5k.»i)  (15) 

Since  Is  a modulo-27r  version  of  we  may  simply  reflect 

all  of  the  conditional  probability  mass  Into  to  obtain  the 

result 

00 

n=-oo 

The  importance  of  this  result  is  the  following:  because  of  the 
Markov  property  of  we  may  write  the  Joint  density  of  (13)  with  the 

conditional  density  evaluated  as  In  (16).  And  as  we  will  find 
in  Section  VI,  (16)  may  be  used  to  assign  transition  probability 
mass  to  a finite  grid  of  phase  values.  These  transition  probabilities 
then  enter  into  the  path  metric  in  a Viterbi  algorithm  for  tracking 
modulo-2ir  phase. 

There  is  an  important  symmetry  property  of  (16):  the  transition 
probability  depends  only  on  the  squared  distance  between  and 
Thus  If  (16)  Is  evaluated  on  an  LxL  grid  of  phase  pairs,  there  Is 
Toeplltz  symmetry  in  the  corresponding  matrix  of  transition  probabili- 
ties. This  means  only  an  L-vector  of  transition  probabilities  need  be 
stored  for  cyclic  reading  to  obtain  the  entries  In  the  matrix. 
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VI.  The  Vlterbl  Tracker  for  Modulo-2n  Phase 
Using  the  result  of  (13),  write  out  the  maximization  problem  of 
(11)  as  follows: 
max  r 


■•k  - * ri  Sc 

k=l  2a  k=l 

Note  that  F,  satisfies  the  recursion 
k 

•■k  ■ ’’k-i  * 2<=k <V*k> 

- llnf(^^)  + 2C^  cos(i|/j^-*^) 

2a 

Thus  the  so-called  path  metric  is 

inff^k/Vi)  + -h  2Ck 

2a 

The  maximization  problem  may  now  be  written 

max  [max  F , + inf(ij,/^„  ,)  + 2C  cos(\|i  -1^)] 

I1 1*"  a 


This  form  leads  to  the  following  observation:  the  maximizing  trajectory 

(call  it  passing  through  on  its  way  to  must  arrive  at 

— — K— 2 

(|i„  , along  a route  {(ji,  },  that  maximizes  F^  - . If  it  did  not  we 
K— 1 K.  1 N— i 

— - - K-2 

could  retain  and  <()j^  and  replace  with  a different  sequence 

to  get  a larger  value  for  F . It  is  this  observation  which  forms  the 
basis  of  forward  dynamic  problem.  The  Viterbi  algorithm,  a forward 
dynamic  programming  algorithm  applied  to  MAP  sequence  estimation  on 
finite-state  sequences,  may  be  applied. 

The  trellis  of  Fig.  1 illustrates  how  the  minimization  of  (20) 
proceeds.  Tabulated  values  of  stored  in  square  array 
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(or  vector  which  Is  cyclically  read)  whose  dimensions  depend  upon  how 
finely  the  Interval  [-ir.n)  Is  discretized.  Call  - x - , with  - ■ 
the  flnlte-dlraenslonal  grid  for  which  defined.  That  Is, 

Is  assumed  to  take  on  only  the  values  1“1,2,...,L,  for  each 

k.  Call  r,  (§,,§  ) the  value  of  the  metric  T,  corresponding  to  a phase 

K.  X m K 

_ Ic 

trajectory  {if,},  which  terminates  at  phase-state  §.  at  stage  k,  after 
J 1 t 

passing  through  stage  at  stage  k-1. 

The  algorithm  begins  with  a computation  of  )^“1»2,...,L 

based  on  measured  values  of  C^  and  \j>^.  See  (18).  If  all  phase  values 

are  equally  likely  a priori,  then  f.nf($j^)  Is  constant  on  all  values 

Otherwise  there  Is  some  a priori  weighting  In  favor  of  some  of  the 

A new  measurement  pair  (C-.ij)-)  Is  obtained  at  k“2  and  T (§  ,§  ) Is 

2.  L z 1 m 

computed  for  m*l,2,...,L  using  a table  look-up  (e.g.  In  ROM)  for  the 

).  The  maximum,  value  of  r-(§,  ,§  ) Is  determined  (the  maxl- 
2 11m  Z 1 m 

mlzatlon  Is  over  all  originating  states  § ) and  the  corresponding  sequence 

m 

A 

(§l^§l)  Is  saved  as  a survivov  sequence  terminating  at  at  stage  2; 

A 

denotes  the  originating  state.  The  survivor  sequence  Is  labelled 
with  Its  corresponding  length  r2(§j^,i2^).  This  calculation  Is  repeated 
for  each  possible  value  of  phase  until  all  pairs  (§^^^,5^2^),  and  correspond- 

A 

Ing  lengths  r2(§£»®jj^)»  1“1,2,...,L,  have  been  computed  and  stored  In 
soft  memory  (e.g.  RAM).  There  Is  a unique  survivor  sequence  corresponding 

A 

to  each  state  £“1,2,...,L.  Caution:  in  the  pair  (§j^,§j^)  the  orlgl- 

A A A 

natlng  state  depends  on  l.e.,  measurements  C2 

and  ^2  discarded  along  with  all  extinct  sequences.  Now  a new 

measurement  pair  (C^.ij/^)  Is  obtained  and  is  computed  for 

mp‘l,2,...,L  using  the  recursion  for  and  a table  look-up  for 


J - - 

^ minimizing  value  of  T2(®j*®j^)«  call  It  r2(§j^,ij^) 
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Is  computed  and  the  corresponding  triple  Is  stored  as  a 

survivor  sequence  originating  In  (*)  at  stage  1,  passing  through  at 
stage  2,  and  terminating  In  at  stage  3.  This  procedure  Is  repeated 
until  L survivor  sequences  of  the  form  determined  and 

stored  with  their  respective  survivor  metrics  • The  measure- 
ments and  all  extinct  sequences  may  be  discarded. 

A A A 

Call  (k)  , ^2  * • • • » ) the  MAP  sequence  estimate 

based  on  k measurements  up  to  stage  k.  This  sequence  has  the  maximum 
value  of  The  parenthetical  notation  (k)  denotes  dependence  on 

measurement  Interval.  In  general  the  MAP  sequence  estimate,  (^^(k+1), 

* . . ’^k+l(^^l)> * measurements  up  to  stage  k+1,  may  differ  from 

the  previous  sequence  estimate  at  every  stage  from  1 to  k.  However,  as 
a practical  matter,  one  can  choose  a sufficiently  large  depth  parameter 
kg  so  that  the  sequence  of  fixed-lag  estimates 

(k)  , k “ k^,  kg+l,...  (21) 

gives  an  approximate  MAP  sequence  estimate.  Here  . (k)  Is  simply 

X-Kq 

the  phase  value,  k^  stages  back.  In  the  MAP  sequence  estimate  based  on 
k measurements.  In  this  way  one  obtains  a phase  track  with  delay  k^. 

Following  Forney  [4]  we  may  summarize  the  storage  and  computational 
requirements  for  the  phase  tracking  algorithm  as  follows: 


(time  Index) 


(®f  » • • •)  . . . ,L  (survivor  phase  sequence 

’ terminating  In  at 

stage  k) 

rj^(i^,i  ...),  1-1,2, ...  ,L  (survivor  metric) 


1-1, . . . ,L 
m— 1 , . . . Ij 


(transition  probability 
matrix) 
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Initialization 

k-1 


“ 5j^  > 1^1,.  ..jL 


ri(§j^,§^)  - lnf($j^-§P  +-^  2C^  co8(*j^-i^),)l-l,2,...,L 

2o 


Recursion 


l « \ 9 * 9 

k m m 


k+l  'I'^k  m' 


•‘"/k+l'SfV 

m 


+ 7T  2Ck+i 

20 

1=1,2,... ,L 


1-1,2,. ..,L 


Measurement/Computation 

C|^  envelope  statistic 

i|)j^  phase  statistic 

Cj^  cos(§j^-i|;^)  + metric 

In  Fig.  1 typical  trajectories  for  this  algorithm  are  Illustrated. 

The  heavy  lines  denote  survivors  and  the  light  lines  denote  path  metric 

calculations  that  are  made  and  then  discarded  In  favor  of  survivors. 

At  stage  3 all  calculations  r.(§„,S  ) are  Illustrated  with  light  lines; 

the  heavy  line  from  to  Illustrates  that  this  path  gives  maximum 

I'o(So>S  ) and  Is  therefore  labelled  a survivor.  Of  course  . 

The  x's  on  the  trellis  Illustrate  sequences  that  have  survived  for  a 

while  before  being  exterminated  by  the  weight  of  evidence.  The  very 

heavy  line  at  each  measurement  stage  k denotes  the  current  MAP  sequence. 

The  circled  number  denotes  the  current  MAP  sequence.  The  sequence  of 

end  points  labelled  with  the  circled  numbers  Is  a sequence  of 

» 

phase  estimates.  Note  this  sequence  of  phase  estimates  differs 
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from  the  MAP  phase  sequence  estimate.  The  latter,  being  a smoothing 
solution,  is  in  fact,  generally  smoother  than  the  former. 

VII.  Simulations 

In  the  simulations  described  here  the  modulo-2n  phase  space  [-ir,iT) 
has  been  discretized  to  M points  ■ J!,2Tr/L-(L-l)Tr/L,  i»0,  1,...  (L-1). 

Figs.  2 and  3,  L=ll.  In  Fig.  4,  L»17.  The  transition  probabilities 
between  these  phase  points  have  been  computed  by  approximating  (16)  on 
an  LxL  grid.  These  values  have  been  stored  in  an  LxL  matrix  (really  an 
L-vector  which  is  cyclically  permuted  to  generate  the  rows  of  the  Toeplitz 
transition  matrix)  and  used  to  calculate  path  metrics  according  to  (19). 

A random  walk  phase  sequeine  was  generated  from  (3)  and  used 

to  construct  the  sequence  of  measurements  given  in  (1).  These  measure- 
ments were  then  used  in  the  dynamic  programming  algorithm  of  Section  VI 
to  generate  the  results  of  Figs.  2-4. 

Fig.  2 is  a relatively  high  signal-noise  ratio  example  where 

a * 0.1  (S/N  = lOdB).  The  phase  variance  parameter  is  a “ (ir/10/T)  , 
n w 

corresponding  to  a phase  trajectory  whose  "1-a  transitions"  are  on  the 
order  of  it/14.  In  Fig.  2 the  solid  curve  with  piecewise  constant  slopus 
is  the  actual  phase  sequence  '^he  dotted  curve  with  piecewise 

constant  slopes  is  the  measurement  sequence  with  ■ arg 

c (-11,11).  The  heavy  lines  labelled  with  numbers  from  1 to  37  are  MAP 
phase  sequences  generated  with  the  modulo-2ii  Viterbl  algorithm  of 
Section  VI.  The  MAP  sequence  after  k measurements  is  labelled  with  the 
corresponding  k value  at  the  end  of  the  sequence.  The  sample  rms  phase 
error  is  approximately  ir/lO,  which  is  about  twice  the  rms  quantization 


In 
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noise  level  of  (ir/ll/J).  Note  In  Fig.  2 that  the  maximum  number  of  "phase 
corrections  to  a path"  Is  4,  occurring  at  measurement  9.  There  is  a 
correction  of  3 phases  at  measurement  27,  in  a region  of  noisy  measurements. 
Thus  for  this  simulation  one  could  have  decoded  the  MAP  sequence  corre- 
sponding to  37  measurements  by  decoding  with  a fixed  delay  (or  depth 

constant)  of  ^^=4.  See  the  discussion  surrounding  (21). 

2 

Fig.  3 shows  a simulation  for  a * 0.5,  corresponding  to  a signal- 

noise  ratio  of -3dB.  The  phase  sequence  is  somewhat  more  wildly 

2 2 

fluctuating  than  in  the  previous  case  : » (n/10)  . The  measurement 

sequence  Is  really  not  underbiased  as  it  appears  to  be.  Noisy  measure- 
ments on  a phase  sequence  that  hugs  the  upper  boundary  of  [-n.n)  are 
bound  to  appeal'  underbiased  because  of  the  way  the  measurement  ■ arg 
e[-Tr,ir)  is  computed.  Several  measurements  In  the  vicinity  of  k * 28-39 
are  reflected  upward  by  2-n  to  further  clarify  this  point. 

The  simulation  of  Fig.  3 illustrates  how  the  modulo-2Tr  Vlterbi 
tracker  acquires  phase  at  low  signal-noise  ratios.  Note  that  early 
on  in  the  run  the  algorithm  jumps  wildly  between  the  phase  points  1 (at 
k=»l),  11  (at  k=2),  9 (at  k=3,4),...,6  (at  k=8) , and  so  on.  Furthermore, 
the  MAP  trajectories  in  this  phase  acquisition  region  are  constant  phase 
trajectories  that  have  high  a pviori  probability;  l.e.,  the  transition 
from  ” §£  to  is  more  probable  than  any  other.  It  should 

also  be  noted  that  the  MAP  sequences  at  k ■ 21-24  are  actually  quite 
good  - the  modulo-2ir  errors  between  the  sequence  values  and  the  actual 
phase  are  small.  Finally,  note  that  for  k > 25  and  26  the  algorithm 
has  been  diverted  from  the  actual  phase  by  a sequence  of  unusually  noisy 
measurements.  By  measurement  k > 27  this  situation  has  been  recognized 
and  a new  trajectory  that  is  constant  through  the  period  of  noisy 
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measurements  has  been  constructed.  This  behavior  is  characteristic  of 
the  algorithm  and  is,  in  fact,  one  of  the  very  nice  features  of  the 
algorithm. 

Fig.  4 illustrates  how  the  algorithm  operates  on  wild  phase  sequences 

that  make  excursions  outside  the  interval  [-x,n).  For  this  example  the 

2 2 

phase  variance  parameter  is  large  : o^  = (if/S)  . In  the  course  of  45 

sampling  instants  the  actual  phase  sequence  wanders  over  a range 

of  3'^  radians.  For  clarity  of  presentation  phase  values  lying  outside 

[-ti,it)  have  been  reflected  Inside  so  that  performance  of  the  modulo-2w  ; 

i 

phase  tracker  may  be  correctly  gauged.  If  loss  of  phase  lock  Is  defined  j 

i 

to  be  a phase  error  of  ir  radians,  then  none  occur.  If  It  Is  defined  to  j 

be  Tr/2,  then  2 occur  in  the  vicinity  of  k=25,  in  a region  of  very  wild  | 

phase  fluctuations.  l 

I 


MIJ IK. 
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VIII.  Applications  to  Data  Communication 
An  appropriate  modification  of  (1)  to  model  data  communication  over 
a perfectly  equalized  channel  Is 


"k  ® ' “k 

Here  may  be  a sequence  of  binary  or  M-ary  complex  source  symbols  which 

may  or  may  not  have  been  trellis  or  convolutional  encoded.  In  the 
simplest  case  tlj^}  Is  an  uncoded  sequence  of  Independent  source  symbols 
drawn  from  an  equlprobable  symbol  alphabet. 

To  model  data  transmission  over  an  Imperfectly  equalized  channel 
with  finite-length  pulse  response  we  proceed  as  follows.  Let  denote 
the  following  (L+l)-tuple  of  phased  symbols: 


+ n. 


(22) 


JVi 

-1  ® . •••  . \_L 


J\-L, 

e ) 


’ \-i 

Denote  by  r'  the  finite  pulse  response  of  *he  channel: 
r'  = (r^,  r^,  ...  . r^^) 


(23) 


(24) 


In  these  definitions  the  superscript  ' denotes  transposition. 

The  noisy,  complex  demodulated  channel  output  at  time  k may  be 
written 

In  the  discussion  below  we  present  two  examples  requiring  succes- 
sively higher  levels  of  complexity  for  simultaneous  modulo-2iT  phase 
tracking  and  symbol  decoding.  The  principle  of  optlnmllty  chosen  In 
these  examples  Is  MAP  data  sequence /modulo-2n  phase  sequence  estimation. 
No  claim  of  minimum  symbol  error  probability  or  minimum  block  error 
probability  Is  made. 


(25) 
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Example  1:  is  a sequence  of  Independent  symbols  drawn  from 

an  M symbol  complex  alphabet  and  transmitted  over  a perfectly  equalized  channel. 
This  example  Includes  such  modulation  schemes  as  M-ary  ASK,  PSK,  QASK,  etc. 

The  number  of  candidate  message  sequences  after  K symboling  periods  is  M . 

2 

The  appropriate  log  density  function  for  this  example  Is 

10, 

- log 

1 K jl  K 

- C - L |z  - I e -h  L log  f(*. /♦.  ,) 

2a  k»l  1^.1  k K X 

It  follows  easily  that  the  path  metric  Is 


(27) 


This  may  be  written  (Ignoring  terms  Independent  of  , I.) 

fC  K 


^ 2Re  z r e (l  P + log  fCL/?.,) 

20  ^ 2a 


An  even  more  explicit  form  of  the  path  metric  Is 


(28) 


^ 2C^  I Ij  co.(*^-i^-e^)  1 |lj 2 + log  f (29) 


where  (Cj^,  i|7^)  are  the  envelope  and  phase  statistics  defined  pre- 
viously and  |lj^|  and  6^  are  the  magnitude  and  phase  (usually  visualized 
in  a complex  signal  constellation)  of  the  complex  symbol  Ij^.  From  all 


three  of  the  forms  for  the  path  metric  It  Is  clear  that  one  is  attempting 

k 


to  reconstruct  phase-corrected  symbols  e 


that  correlate  well 
with  the  received  data  while  also  satisfying  a smoothness  constraint 


(as  embodied  In  the  j^))  mod-27r  phase  sequence. 

2 

This  Is  actually  a mixed  density-probability  mass  function.  However,  be- 
cause the  probability  mass  function  for  assigns  uniform  mass  M~  to 

each  K-tuple  we  Ignore  this  term  and  call  f (•,*,*)  a density. 
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The  appropriate  trellis  for  Viterbi  decoding  of  the  niodulo-2Tr 
phase  sequence  and  the  symbol  sequence  Is  the  MxM  symbol-  phase  trellis 
of  pairs  m=l,2,...,M,  )l=l,2, . . . ,M,  corresponding  to  all  allow- 

able pairs  of  symbol  and  phase.  Recall  M denotes  the  number  of  phase 
points  in  the  discretized  modulo-2'ir  phase  space.  If  is  drawn  from 
a binary  alphabet  for  all  k,  then  the  symbol-phase  trellis  simply  con- 
sists of  2M  symbol-phase  pairs  over  which  survivor  sequences  are  formed 
as  outlined  in  Section  VI. 

Example  2:  ^ sequence  of  independent  symbols  drawn  from  an 

M symbol  complex  alphabet  and  transmitted  over  a channel  with  finite 
pulse  response  For  this  example  the  appropriate  log-density  for 

the  K measurements  is  (neglecting  uninteresting  constants) 

log  f({z^}^.  {I,}^)  = 


(30) 


K K 

2a  k=l  k=l 

Recall  Xj^  depends  on  a specific  set  of  source  S3nnbols  and  phases  as  in 
(27).  If  the  modulo  2Tf  phase  space  [-tt.tt)  is  discretized  to  M values 
then  there  are  M distinct  vectors  Xj^. 

The  path  metric  for  this  problem  is  simply 


2a 

The  trellis  for  Viterbi  decoding  consists  of  the  M allowable  X- 

vectors  corresponding  to  M phase  values  and  distinct  symbol  strings 

of  length  L+1.  The  path  metric  for  each  allowable  transition  in  this 


trellis  is  computed  as  in  (31)  and  used  to  construct  survivor  sequences. 
Of  course,  not  all  transitions  are  allowable  because  specification  of 
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(I  T T ^ restricts  the  next  allowable  sequence  to  one  of 

M sequences  beginning  with  terminating  in  through 

This  means  one  must  store  H survivors  (“M  2^^^  for  binary  signalling) 

L+1  ^ L+2 

and  compute  path  metrics  for  M M M.  M allowable  transitions  (“M  2 

for  binary  signalling) . 


IX.  Conclusions 

We  have  posed  the  problem  of  obtaining  a MAP  sequence  estimate  for 
phase,  modulo-2ir.  To  solve  the  problem  we  have  derived  an  expression 
for  transition  probabilities  in  a modulo-2'iT  version  of  random  walk  and 
used  the  result  to  derive  a dynamic  programming  algorithm  for  MAP  sequence 
estimation.  Simulations  indicate  the  performance  of  the  algorithm  la 
excellent.  Computational  demands  are  modest,  making  the  appraoch  attrac- 
tive for  data  communication  applications. 
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